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1. INTRODUCTION 


The present paper extends a previous analysis [1] of weak rotating 
waves in a centrifugal compression system with a vane less diffuser, shown in 
Fig. 1.1 adapted from that paper. Conditions were found for the existence of 
permanent rotating waves of small amplitude. The theory assumed an 
inviscid flow process in the vaneless diffuser. Waves of "fast" and "slow" 
tvpes were found, consistent with experimental experience. 


The slow waves (usually called "diffuser stall" in the literature) are 
especially interesting in that resonance was found. That is, when flow 
coefficient <t> is reduced, a value is suddenly reached for which pressure 
amplitude may be arbitrarily large, for a given radial-flow disturbance of a 
given mode. Fig. 1.2 from [1] shows the real part of a presumed pressure 
disturbance, and the corresponding slow-wave speed coefficient. Negative 
values of the indicated function presumably represent amplification, but that 
could not be proved in the linearized analysis. This predicted resonant 
behavior is consistent with experiments, and is quite unlike anything found in 
the corresponding problem of axial system stall [2]. The present paper 
begins with remarks showing how the results of [1] can be reconciled with 
experimental evidence [3] of an effect of diffuser width on the occurrence of 
rotating stall. 

The weak-wave analysis cannot, of course, predict limit cycles, nor 
can it predict what will actually happen when flow coefficient is reduced 
through the resonant condition. The wave cannot remain weak and 
permanent (the only permitted linear solution); neither can its amplitude be 
infinite. It must grow, and yet disturbances must physically remain finite. 
Clearly, the next step should be a nonlinear analysis, and the appropriate 
approach would seem to be that of [4] . 

Therefore, the pressure and mass balance equations for the centrifugal 
system are first developed, as in [4]. A Galerkin procedure is applied, also as 
in [4], for an assumed single rotating mode, and the corresponding dynamic 
pressure rise in the diffuser is found to second order in disturbance 
amplitude. Subject to certain assumptions to be described, a final set of 
equations is derived, which govern a rather general transient comprised of 
surge-like and rotating-wave disturbances in combination. In the solution of 
these equations, use may be made of the linear solutions for velocity and 
pressure disturbance already found in [1]. 

Then, an ideal example problem is solved with these equations, and it is 
shown how resonance is traversed when flow coefficient is decreased 
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(throttled). In effect, a surge transient must generally occur near resonance, 
so that the wave amplitudes, though singular, remain finite. 

Away from resonant conditions, the possibility of more than one wave 
type should be considered, because two families of rotating waves (fast and 
slow) are in principle allowed. This problem is discussed. 
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2. EFFECT OF DIFFUSER WIDTH 

The assumption of inviscid flow in [1] and this Report may be 
questioned in view of the known effects of Reynolds number on performance 
and stability of the machines under discussion [3,4]. In most flow problems, 
however, it is important to know what inviscid analysis has to tell, that 
knowledge avoids the pitfall of ascribing puzzling flow observations, 
especially instabilities, to viscous phenomena such as separation, when the 
real cause may lie elsewhere. It is believed that the present study shows an 
inviscid basis for rotating instabilities found with vaneless diffusers. 

Diffuser width has no effect in the present theory, yet experiments, 
reviewed and analyzed in [3], show a clear tendency for increased diffuser 
width to increase the critical flow coefficient for instability. Fig. 35 of [3] 
indicates that when width / increases from 2 percent to 15 percent of exit 
radius if, the critical flow coefficient increases from about 0.10 to about 
0.25. An asymptote of perhaps 0.28 is apparently approached as / is 
increased further. Fig. 1.2 (Fig. 10 of [1]), on the other hand, gives a critical 
flow coefficient of about 0.3; thus the inviscid analysis agrees with the 
asymptotic value for large l suggested by [3]. The influence of radius ratio 
will be different in [1] and [3], because the discharge is different. In any case, 
a larger radius ratio leads to a higher critical flow coefficient in both theory 
and experiment. The variation with width must obviously be an effect of 
viscosity, which was not considered in [1]. 

Fig. 36 of [3] also suggests that the foregoing critical flow coefficients 
are valid for Reynolds number (based on / and on impeller exit velocity) of 
10 5 , but should be increased by about 0.07 for 10 4 and decreased a like 
amount for 10 6 . 

Boundary-Layer Displacement and Diffuser Convergence 

To explain these results in terms of the present theory, one should first 
consider the displacement effect of the wall boundary layers in the diffuser. 
The sketch (Fig. 2.1) shows a narrow and a wide diffuser, each with the same 
growing boundary layers and hence the same displacement effects. Clearly, 
the narrow diffuser will effectively be more convergent because of the 
boundary-layer displacement effect. 

Convergence effect was discussed in [1], and was expressed as a power 
law (See Fig. 9 of [1]) 
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Results shown in Fig. 10 of [1] showed that the critical flow coefficient 
(in the first circumferential mode) was lower for a convergent case (n= l) 
than for the parallel-wall case (n = 0). Subsequent calculations give the same 
result for the higher modes (3rd and 4th) which are the most unstable. Figs. 
2.2 and 2.3 show the pertinent results. Thus, convergence due to boundary- 
layer displacement would be predicted to cause a lowered critical flow 
coefficient for a narrow diffuser compared with a wide one. Of course, 
when the width is much greater than the boundary-layer thickness, further 
increasing the width will have no effect. One may consider that the present 
inviscid theory assumes infinite diffuser width, in effect. 

The foregoing discussion is qualitative; to make numerical predictions, 
or to explain the Reynolds number correction also suggested in [3], would 
require quantitative knowledge of the displacement thickness and how it 
varies with Reynolds number. That information cannot readily be inferred 
from the literature. 

The theory predicts that a more convergent diffuser will be more 
stable, whether the convergence is due to displacement effect or is simply a 
geometrical feature of the diffuser. It would therefore appear advantageous 
to build the diffuser with a convergent shape. However, pressure rise in such 
a diffuser would be more gradual, and to achieve the same final pressure 
ratio, its exit radius would have to be larger. A larger exit radius, in itself, is 
destabilizing. For the third circumferential mode and parallel diffuser walls 
(n = 0), the critical flow coefficient is 0.180, if rF is 1.6. If the walls are to be 
convergent to the degree that n = 1 , but the pressure-rise coefficient is kept 
the same, then one finds that ip must be 1.89 instead of 1.6. For tp = 1.89 
and n= 1, the critical flow coefficient turns out to be 0.175, which is virtually 
the same as that for the shorter, parallel-wall diffuser. Therefore, it would 
seem that if pressure coefficient is held constant by changing exit radius, 
convergence of the diffuser has negligible effect on stability. 

Details of Linearized Analysis 

Owing to length restrictions, equations and formulas which included 
the effects of convergence and mode number were not presented in [1]. For 
the record, Appendix A of this Report collects those details, thereby 
generalizing the Eqs. (9 through 20) which appear in [1]. The relevant 
discussions and arguments appear in [1] and need not be repeated in the 
Appendix. 
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3. TRANSIENT EQUATIONS 

Pressure and Mass Balances of System 

The overall pressure balance of an axial compression system appears 
in Eq. (42) of [4], and is repeated here for reference: 

¥© + / C <D^ T c (<D + g) - mY^ + ^ (2Y^ee + Yee) (3.1 ) 

For the centrifugal system of present interest, shown in Fig. 1.1, the 
corresponding equation may easily be derived: 

y + i c = 'PcCO+g) - mY^ + ^ (Y^ee + Yeee) - [dPl + (Pi - Pgsl)] (3-2) 

(a) (b) (c) (d) (e) (f) 

The terms of this equation are lettered, and will be discussed in turn, making 
reference to the corresponding terms of Eq. (3.1). 

(a) The total-to-static pressure-rise coefficient of the entire system is 

¥ = (PF - PT)/pU 2 (3,3) 

where pj is considered constant, while the plenum pressure pF to which the 
diffuser discharges may be a function of time. The rotor tip velocity is 
denoted by U. 

(b) As in Eq. (3.1), the second term describes flow inertia up to the rotor 
exit, based on the coefficient d> of circumferentially-averaged radial velocity 
there. The dimensionless length l c is now 


, , Aft J_ 

/c-/l Afo + 2a’ 


a = ri/xU 


(3.4) 


Unlike the axial case, the inlet velocity differs from the rotor exit velocity 
used as reference, and, by continuity, goes inversely with flow area Af. 
Therefore, the area ratio Afi/Afo now appears. If x is defined as in Eq. (2) of 
[4], then a new factor of 1/2 appears in Eq. (3.4) because there is no stator in 
the present case, whereas [4] assumed a two-row (rotor-stator) stage. The 
time coordinate is relative to time for the rotor to travel one radian. 
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(c) As usual [4], the axisymmetric, quasi-steady pressure-rise coefficient is 
denoted by 4A, stipulated to include the pressure rise in the diffuser, which is 

written as - Pqsl) 

(d) As discussed in [4], the function Y(^,0) is the potential for any wave-like 
velocity disturbance, depending only on ^ and 0; such a disturbance is 
assumed "straight through", not a function of position along the flow path 
from entrance to exit. Corresponding to [4], the disturbance of radial flow 
coefficient at exit is denoted by g, and 

(|) = 0 + g ; g = -Yee (3.5) 


The term (d) specifically refers to the pressure change caused at the entrance 
by the appearance of a transverse (that is, circumferential) velocity 
disturbance (w) there, which is in turn induced by the axial disturbance 
described by g. The assumption of [4] is carried over, that w and g are 
related as if the flow were potential and harmonic at the entrance. The 
applicable relations are 


Aft 

w 1 = U XS Y 9 ; 


m = 


X An 

tF Afo 


(3.6) 


Again, the area ratio Af i/Afo appears, because the disturbance function g is 
defined at the rotor exit, while Eq. (3.6) refers to events at the entrance. The 
radius ratio ip also appears as a scaling for the potential function at the 
entrance, because all distances have been referred to the exit radius. 

(e) This term is the same as the corresponding term of Eq. (3.1) with a factor 
of 2 omitted. That difference reflects the presence of only one blade row, 
while the axial case [4] has two. 

(f) The last group of terms is new to the centrifugal case, and refers to the 
pressure in the vaneless diffuser. A pressure coefficient in the diffuser is 
defined below, and may be assumed to consist of an angle-averaged part 

(perhaps unsteady) P and a wave-like part dP (whose angle average is zero). 


plj 2 = P = P(r,4) + dP(r4,0> 


(3.7) 


6 


• ii in in 



Now, the axisymmetric overall pressure rise function \]/c , presumed known, 

includes quasi-steady pressure rise - P qsi (d>) from diffuser entrance to exit. 
Hie re fore that quantity must appear in the group (0, so that Eq. (3.2) is 
satisfied in the axisymmetric quasi-steady case, when 4* = V c and all its other 
terms must vanish by definition. 

In the postulated system (Fig. 1.1), pressure falls from plenum to 
exhaust through a throttle. The throttle has a characteristic performance 
function Ft, presumed known, just as does the compressor: 

v F = Ft(Ot) (3.8) 

The flow coefficient Or may differ from O because mass may accumulate in 
the plenum. In fact, with the plenum regarded as an isentropic air spring, 

Eq. (36) of [24] applies: 

/c'f'E = ~ ! T( < I>- < t>r) 0.9) 

4Bg 

Eqs. (3.2), (3.8), and (3.9) constitute the governing equations for 
disturbances to be considered. The problem is to find permitted solutions for 
<£>(%), Y(%,0), and corresponding pressure distributions in the diffuser. 
First, of course, the throttle and compressor characteristics must be 
specified. 

Galerkin Version with Cubic Characteristic 

As in [4], the system is greatly simplified by use of the Galerkin 
procedure, which represents the periodic functions Y(£,0), and for the 
present problem, also dP(r,^>9), as harmonic functions with only their time- 
dependent amplitudes remaining unknown. It is assumed that Y(or the 
related g) can be well represented by waves traveling at dimensionless speed f 
in the positive 0-direction; thus, as in Eq. (58) of [4], 

g = W Ai cos 0 (3.10) 


and 
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P = P (r,£) + jP (r,^) cos 0 + 2 P(r,^) sin 0 


(3.11) 


where two new functions of radius and time, iP and 2 P have been used to 
construct dP, and where 


0 = 0-f£ (3.12) 

Two major questions may be raised about this use of the Galerkin 
method. First, there may be more than one type of wave represented at a 
given time; Ref. [1] shows that there are "fast" and "slow" wave types and 
also higher harmonics in the centrifugal, but not the axial, case. These 
complications can easily be dealt with in the present analysis, as will be shown 
later. The second question is more fundamental; waves in this nonlinear 
problem may differ substantially in shape from the simple harmonic form. 
Fig. 5 of [4] shows this difference for the axial problem. No such 
comparison can be made for the centrifugal case; the necessary exact 
nonlinear solution of the diffuser flow is not available. However, wave 
traces measured in centrifugal machines [3] seem to look rather more 
harmonic than is usual in axial machines, perhaps because of a smoothing 
effect of the diffuser. Therefore, one may be optimistic about the accuracy of 
the Galerkin method for the centrifugal problem. 

The Galerkin method is especially needed in the present problem 
because of the need to provide a nonlinear solution for flow in the diffuser, 
whereas in the corresponding axial problem, the only nonlinearities arise in 
the compressor and throttle characteristic functions. 

It is further assumed that the compressor characteristic is well 
represented by a cubic function, namely 


Vo = ¥o + H[| ^-1(^)3] 


(3.13) 


This function, shown in Fig. 3.1, has one more parameter (S) than the 
corresponding cubic for axial systems (Fig. 3 of [4]). Because measured 
characteristics for centrifugal machines tend to be very smooth functions, the 
assumed cubic representation should be very satisfactory. The throttle 
characteristic might be a simple quadratic, also illustrated in Fig. 3.1. 

To carry out the Galerkin procedure, one first substitutes Eqs. (3.5), 
(3.10), (3.1 1), and (3.15) into Eq. (3.2), and then performs three 
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integrations over a cycle of 0, namely of Eq. (3.2) itself, of Eq. (3.2) 
multiplied through by cos 0, and then of Eq. (3.2) multiplied through by 
sin 0. The results of these operations are, respectively, 

4 / + / c Op = \|/ C (^>)-|^( < J > -S) Ai - (Pi -Pqsl) ( 314 ^ 

(m4)A, 5 + | {§A?-f[l -<^]A,} (3-15) 

[( m + ^) f -^]Al+^2 p l© = ° 0.16) 

These equations are nearly equivalent to Eqs. (59 through 61) of [4], The 
next step must be to analyze the diffuser flow in order to relate the new 
pressure terms to velocity-disturbance amplitude. 

Pressure Disturbance in Diffuser 

Flow in the vaneless, parallel-wall diffuser will be considered inviscid, 
just as it was in [1]. The governing Euler equations, written in original 


dimensional variables, are 


3 1 n u/2 WVo 

Vl+ 3r [;2 (v2 + w2) + p 1 -(T + WWr " r ) = 0 

(3.17) 

Wt + r ^[2 (v2 + w2) + p ] + (7~ +VWr ' r ) = 0 

(3.18) 

v 1 n 

Vr + 7 + 7 W0 = O 

(3.19) 


Next, the Galerkin definition of radial velocity disturbance, Eq. 

(3.10), is generalized to allow for phase change in the diffuser: 

v = U + WA(r, £,) cos 0 + WB(r, ^) sin 0] (3.20) 
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Using this equation in the continuity equation, Eq. (3.19), one may integrate 
the latter to find circumferential velocity: 

w = U [^ - W (A + rA r ) sin 0 + W(B + rB r ) cos 0 + X(r£) (3.21) 

The first term implies rotor blades are radial at rotor exit. The function X 
represents a swirl of unknown intensity, resulting from integration. The 
corresponding pressure representation is already in place, as given by Eq. 
(3.11), with Eq. (3.7) providing the connection with static pressure. The 
pressure in the diffuser, Pp, may be a function of time, of course, but will 
disappear when differentiated in Eqs. (3.17) and (3.18). 

The applicable boundary conditions in the diffuser are as in [1]. Radial 
velocity is specified at the difftiser entrance: 

A(l£)sAi©; Bi=0 (3.22) 

Circumferential velocity disturbance vanishes there, and therefore 

Xi=0; A ri = -Ai; B r , = 0 (3.23) 

At the diffuser discharge, the pressure disturbance must vanish: 


p F = iPp = 2 Pf = 0 (3.24) 

Two further assumptions are now made about the diffuser flow, in 
order to facilitate solutions; first, that propagation speed f is constant (In 
problems involving inlet distortion, f is no doubt a function of time), and 
second, that terms of order the cube of the given disturbance A may be 
neglected; this means that the square of the unknown swirl X, and products 
such as AX, may be neglected. The equations resulting from the Galerkin 
substitutions into Eqs. (3.17) and (3.18) and integration after multiplying by 

1 , cos 0 , and sin 0 are six in number: 

Pr = -i<p^ + (l + <t>2)i + i r W2(Ar + B^) + 2W2y (3.25) 


lP 2 = W { -A ? - ® (A r - 7 ) + [ 5[2rB r + (1 + f r 2)B] } (3.26) 
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(3.27) 


2 p r = W { -B^ - J (Br - 7 ) - jz U rA r + ^ + fr2)A J 1 

0 = 3(AB r - BA r ) + r (A Brr - BArr) + 2 [X^ + ^ ( rX )r] (3.28) 

= W [-r(A + rA r )^ - (1 - fr 1 2 ) (B r + ~) - O (rArr + 3A r + “)] (3.29) 

2 P = W [- r (B + rB r )^ + (1 - fr 2 ) (A r + ~) - <*> (rBrr + 3B r + “)] (3.30) 

Even with the Galerkin simplification, the foregoing system involves 6 
unknown functions (3 for pressure and 3 for velocity) as well as an 
eigenvalue (f) to be found. 

The desired solution is permitted by a final assumption, namely that 
the diffuser flow is quasi-steady. That is, one assumes that surge-like changes 
in the overall system will have a much longer characteristic time than the 
time required for fluid to traverse the diffuser. For this to be true, /c must be 
much greater than 1 . There may be cases when this is not true, of course. 
Under this assumption, all derivatives with respect to £ may be dropped, and 
three benefits become apparent: 

1 .inear wave . Eqs. (3.26), (3.27), (3.29), and (3.30) are now identical 
to those which govern the linear, small-disturbance problem of [1]. 

Therefore, properly interpreted, the pressure-wave and velocity- wave 
solutions of [1], including resonant behavior, can be carried over directly to 
the present nonlinear problem, Eqs. (3.14 through 3.16). Nonlinearity 
resides only in Eqs. (3.25) and (3.28), which now can be used separately to 

find P and X. 

Swirl -function solution . Eq. (3.28) can be integrated to find the swirl 
function X . The result is 

X = ^r2(BA r -AB r ) < 3 - 31) 

A constant of integration vanishes in view of applicable boundary 
conditions, Eqs. (3.22) and (3.23). This result is consistent with the previous 
supposition that X is second order in velocity amplitude. 
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Mean pressure-rise solution . Knowing A, B, and X, Eq. (3.25) can be 

integrated to find the mean pressure rise function needed for Eq. (3.14). 
One finds 


P = 




[r(AJ: + B 2 ) + |(BA r -AB r )]dr (3.32) 

r F r 


Integration limits are chosen to make Pp vanish, as Eq. (3.24) requires. 
Obviously, in the quasi-steady, axisymmetric limit Eq. (3.32) reduces to 

Pqs = - 2 ( 1 + < t> 2 ) ( 72 " Z~2* (3.33) 

c r r p ^ 

For use in the last term of Eq. (3.14), one may write 


1 r F y 

Pi - Pqsl = - y W2 J [r (A? + B r 2 ) + s (BA r . AB r )]dr 


<D 


(3.34) 


Final Equations for Transient 

To write the final transient equations, it is necessary to interpret A , B, 
lP, and 2?* defined in Eqs. (3.1 1) and (3.20), in terms of quantities found in 
the small-disturbance analysis of [1]. In that treatment, radial velocity was 
defined in complex form 

7j = ^ + Ae$; A = AR + i AM (3.35) 

and, for convenience, 

ARi = 1 ; AMi =0 (3.36) 

Therefore, in terms of the known linear solution, Eq. (3.20) should be 
written 
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^j = Y + WAj (A cos 0 + 6 sin 0) 

with the understanding that 

A = AR ; B = - AM (3.38) 

In [4], and also in the present analysis, it will prove convenient to express Ai 
in terms of its square: 

J = A? (3.39) 

In [1] the complex pressure rise coefficient had the form 

dP = (PR + iPM) e*® (3.40) 

and therefore, in terms of the previous linear results, Eq. (3.1 1) should be 
written 

d P = WAi(PRcos0-PMsin0); PR = PR; PM = PM (3.41) 

Therefore, where one sees A, B, iP, and 2 P in Eqs (3.14 through 

3.16), one should write instead WVT A, \WT^, WJ ffo, and -WVTPM . 
This done, one may write the final transient equations as follows: 

+ l c <X>^ = y c (O) + 

|{W2 | F [r (A r 2 + ^) + |(M r A6 r )]dr-|^(<l>-S)) J (3.42) 


(m 


+ ^» J ^ + wi4 J - 3 h-( ; 


W 


W A . 

)2] + 2^-PRi} J = 0 



W/H r 0) _ 
2 L\V 

4B g 


FjC'F) 

W 



(3.43) 


(3.44) 
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These equations parallel Eqs. (59-61) of [4], but differ in the way wave 
amplitude appears. For example, pressure amplitude affects rotation speed f 
in the present case. Most notably, Eq. (3.43) shows that a negative value of 
the real part of the pressure wave (the part in phase with the radial velocity 
disturbance), in itself, causes amplitude (J) to grow with time. This confirms 
expectations of [1]. As one expects from [4], the same equation, together with 
Eq. (3.13), shows that positive characteristic slope also makes wave 
amplitude grow. 

The most interesting question for these equations to answer is that 
posed in the INTRODUCTION, namely, how the known resonance 
properties of the linear solution come into play during nonlinear transients. 
The section to follow will focus on this issue. 
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4. TRANSIENT SOLUTIONS 


Singularities of Linear Solution 

A A A 

Ref. [1J showed how A, B, and PR have singularities as functions of 
How coefficient O. It can easily be shown that these singularities are all 
simple poles in general. For the present purpose of illustration, a particular 
case will be considered that is easily analyzed, namely parallel walls (n = 0), 
rp = 2.08, and the fundamental mode (b = 1), for which Fig. 1.2 suggests that 
propagation speed f vanishes at the resonant condition, <t>* = 0.15464. 
(Asterisks will denote the resonant, or singular condition.) As flow 
coefficient is changed, the value of f will also change, but as Fig. 1 .2 
illustrates, f is minimum at resonance. Thus it is surely accurate to assume f 
constant when <!> is near O*. 


For the particular case defined above, 

f = 0 ; n = 0; b=l (4.1) 

the solution of Eqs. (9 through 1 1) of [1] can be written, in notation of [1], 


A = ar 1 ' l/<I> + Pr - 2 + y 

(4.2) 

w = i [oc(2 - r 1 ' ^ - |3r -2 + y] 

(4.3) 

rP = -2<Da(2 - -j~) rL + (i - <D) 0r -2 - (i + O) y 

(4.4) 

The boundary conditions 


A(l) = 1 ; W(l) = 0 ; P(rp) = 0 

(4.5) 


allow (3 and y to be expressed in terms of a , and provide the real and 
imaginary parts of a as follows: 

G -2 M .2 

areal = Q2+M 2 ^ ( 1+ r F ) + G 2 +M 2 r F ^ 
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G _7 M -2 

aimag = G 2 + M 2 ® r F ) + G 2 +M 2 ^ ( 1+ r F^ 


(4.7) 


where these definitions apply: 

In r F 

* = ( 4 - 8 ) 

l-O 2 1+30 2 

G = - 4rF sind - 8 OrF cosO + — + — :r — (4.9) 

Or F v 

M = 4rF cosd - 8 OrF sind + 2 (1+ rj?) (4.10) 

Eqs. (4.6) and (4.7) indicate that there is a singularity only when G 
and M vanish separately, and solving Eqs. (4.9) and (4.10) gives the 
conditions 


(3 4-1) (D-Z 

cosd* = t (4.11) 

2(l+4<h* 2 )rp 


1 +4+<D* 2 (3+7rp) 

sind* = ~ (4.12) 

4<D*rp(l+40>* 2 ) 


which in turn can be solved to discover that, for this case (Eq. (4.1), 
resonance occurs at 

<£* = 0.15464 if r F = 2.08 (4.13) 

The singularity is clearly a simple pole, and one may use the foregoing 
results in Eq. (4.4) to find the real part of the pressure singularity, needed 
for Eq. (2.43). Because G and M vanish there, derivatives with respect 
to d>* will appear: 
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(4.14) 


G'(0*) = -4 rF [(- + 2) cosi}* - 20*(-^») sintl* 

1 — t[ 1+4>*2+ (1-30*2) fp] 

0* 2 i> 

/"v (V 3|C 

M’ (O*) = -4rF [(- + 2) siird* + 20* (- ^r) cosd*] (4.15) 

In present notation, then, the pressure singularity is 

1*1 = o W t ^'[(1-30*2) (1+ r"p) -4(1- r F 2 )] 

+ ^ [(1-30*2) (i. r[ ?) + 40*2 (1+ r F )] } (4.16) 


It is necessary also to find the singularity of the integral in Eq. (2.42): 


1= [r(A r 2*B r 2) + §(M r -M r )]dr (4.17) 

1 

After arduous but straightforward analysis, one finds the result: 


1 

1 “ ( 0 - 0*)2 


0* 2 (1+ r’p) 2 + (1- r’p) 2 r 
[ L’ 2 + M' 2 ^ l f‘2 r F 2( 0^2' l) + 


\ ( 1+ r p) ’ JJa- r F) + [ r p ' 0* 2 (rp ' r F ^ COsd * 

-2 - (3nF + ~) sintf* } (4.18) 

This quantity, it should be noted, is a second-order pole. 

Inserting Eq. (4.13) into Eqs. (4.16) and (4.18) gives numerical values 
for pole strengths in this particular case: (For simplicity, 0 will hereafter be 
measured from the appropriate critical value; that is, 0* = 0.) 
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1 

m + ^ 1 

2a 1 


A 

PRl = - 


b* 


2/c 


CD ; 


where b* = 0.020166 


2/c 

1 

+ 2a 


= -2c* 


02 ’ 


where c* = 0.0005062 


(4.19) 


(4.20) 


Simplified Example Near Resonance 

Reduction of equations . It is now possible to illustrate the behavior of 
the transient equations (3.42, 3.43 and 3.44) near the resonant point for ip 
= 2.08 (chosen for the analytic simplicity which results). Only the singular 
parts of the coefficients of O need be kept. Also, because the focus is to be 
on the effect of resonance, the characteristic function \j/ c (<I>) will be assumed 
locally to have a constant slope SL. Zero slope would imply neutral damping 
for disturbances. Arbitrarily, one may then choose \j/ c (0) = 0, and also 
measure flow coefficient from the resonant value. The three transient 
equations then reduce to 

x F-SLO + / c <^ + c*^2 = 0 (4.21) 


U4-b*£ = 0 

(4.22) 

k'Vl- 1 , (<S - <t>T) = o 
4B g 

(4.23) 

These equations can further be simplified 
following scaling transformations: 

in appearance by making the 

0) = 2B g b*q 

(4.24) 

T = b*z 

(4.25) 

c* 

(4.26) 

§ = 2B g ; c C 

(4.27) 
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(4.28) 


SL = a/2B g 

so that the three governing equations become 


A. n 

z-aq + q^ + ^2 = 0 

(4.29) 

o 

ii 

<-<1 cr 

l 

(4.30) 

z^ - (q - qT) = 0 

(4.31) 


At C, = 0, initial values of the quantities z, X, and q (which must be near 0, of 
course) would be specified, and then would be caused or allowed to change, 
perhaps under influence ol the throttle coefficient qT- 

Small perturbations . Because q must be and remain small for the 
governing equations to be valid, it seems reasonable to further transform 
variables to reflect a small initial value of flow-coefficient disturbance, e: 


q = eQ ; z = eZ ; ?i = e 2 L (4.32) 

The governing equations now are 

8Z-eoQ + eQ^ + ^j2 = 0 (4.33) 

eL^_^ = 0 (4.34) 

Z£ - (Q - Qt) = 0 (4-35) 

with initial conditions 

Q(0) = 1 ; L(0) = L 0 ; Z(0) = 0 ; Qt(0) = Qt 0 ( 4 -36) 


It is necessary to specify an initial rotating-wave disturbance level 
(Lo), which might arise from inlet distortion or some similar cause. 
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Although e may be considered much less than 1 , the equations do not 
thereby become linear, because of the singularities. However, the equations 
may be expanded in powers of e, and "inner" and "outer" solutions may be 
identified. These solutions will be identified and discussed, and then related 
numerical solutions will be presented. 

Inner and outer solutions . The appropriate "inner" variable is 

c = C/e (4-37) 

and when this variable is used, the inner small-e solution proves to be 


o 

II 

N 


(4.38) 

2 + Ifa|S±| 
k + 2 m 1 Q+k 1 

II 

N1 1 1—1 

(4.39) 

L=}(k2-Q2) 


(4.40) 


where k is a constant related to the initial value of L : 

k2=l+2Lo (4.41) 

The flow coefficient Q defined by Eq. (4.39) and the corresponding 
rotating-wave amplitude L defined by Eq. (4.40) are displayed on Fig. 4.1. 

This inner solution provides the most significant result of this paper, 
that flow coefficient cannot decrease smoothly through a resonance value of 
<I>. Rather, an abrupt drop must occur, to a level just as far below resonance 
as it was above resonance when the transient began. During this drop, 
rotating-wave amplitude increases to a finite maximum level 

Lmax = 5k2 = !+Lo (4.42) 

and and then decreases. It is noteworthy that this increase is limited to a value 
of 0.5 above the initial level. 

The beginning of the process just described can be imagined to occur 
as the result of a reduction of throttle flow Oj, which brings O into the range 
of negative damping; above that level, the terms involving L would not 
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appear. Fig. 1 .2 shows that amplification begins, for the case of present 
interest, when O is about 0.01 above resonance. Consulting Eqs. (4.24) and 
(4.32), and assuming that B g is about 1 and l c is about 10, one finds that the 
choice e = 0.01 places the beginning of the transient (Q = 1) at the beginning 
of amplification. That is the value to be used subsequently. It is interesting 
that neither the characteristic slope nor the throttle schedule have any 
influence on the abrupt, "inner" process. 

The influence of initial wave amplitude Lo may be inferred from Eq. 
(4.39) or Fig. 4.1. Especially, the limit as Lo vanishes is of interest. In that 
limit, k approaches 1 , and the initial conditions Q = 1 and Q = 0 will be met 
only at = -<*>, on the "inner" scale. This reasonable result means that if the 
initial wave amplitude is very small, then the time for the wave to grow will 
be very long, but only on the short time scale of the drop to follow. 

One may question why the singular solution goes only one way; why an 
increase of d> is not allowed. Apparently, d> can increase smoothly 
through resonance because it approaches resonance through values of O for 
which the relevant wave is strongly damped. Then as it reaches values above 
resonance, the amplification has no wave amplitude to amplify, and the 
amplitude may remain zero as <f> continues to increase. 

The "inner" time scale is very short; if e = 0.01, shorter by a factor of 
100 than the "outer" scale which is based on £ itself. The outer solution does 
depend on the throttle and the characteristic slope; if one is held constant and 
the other set equal to zero, the leading outer approximation is simply a 
harmonic surge wave with an amplitude matched to the final (negative) inner 
limit for Q, with no rotating-wave amplitude, and with phase determined by 
matching to the inner value of Z at £ =0: 


L = 0 (4.43) 

Q = Qt - (Qt + k) cosC, (4.44) 

Z = - (Qt + k) sin£ (4.45) 

Further consideration of matching shows that, to leading order, the 
outer and inner solutions have the same origin, namely where Q crosses 
resonance. To first order in e, a phase shift would apply in the outer 
solution. 
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Thus, a surge transient must follow the sudden drop at resonance, 
triggered by that drop. The subsequent course of the surge transient 
naturally depends on the throttle schedule and the characteristic slope. In 
principle, however, at any subsequent time when the surge transient carries 
downward through resonance, another sudden change with rotating- wave 
features could occur. 


Numerical Solutions 

Eqs. (4.33 through 4.36) can easily be solved numerically, of course, and a 
number of such solutions for various parameter choices will next be 
described. The foregoing theoretical results show that the "inner", abrupt 
changes will depend only on initial conditions of Q and (slightly) L. The 
throttle schedule and characteristic slope will have effect only on the "outer" 
scale; therefore, the figures will present results on that scale (Q. 

Slow throttle closure . Fig. 4.2 shows results (on the "outer" scale) of a 
numerical solution of Eqs. (4.33 through 4.36), for the indicated parameter 
choices, which correspond closely to the theoretical results just discussed. 
Amplification was assumed to occur only in the range 1< Q < 0, and damping 
only in the range 0 < Q < -4, as Fig. 1.2 implies. The essential action 
indicated on Fig. 4.2 results from a slow throttle closure (Ki) from an initial 
level (Ko), just above the amplification range, where both Q and Qt agree. 

Qr = K 0 -KiC ( 4 - 46 ) 


The results illustrate the "inner and outer" solutions; the fl ow 
coefficient slowly descends, and pressure-rise function Z slowly rises, in 
response to Qt , until the amplification range is entered at Q =1. Then a 
sudden drop (Eq. (4.39)) occurs, together with a burst of rotating-wave 
amplitude (L). Following that drop, a surge process begins and continues 
indefinitely (The characteristic slope a is zero in this example.) 

As anticipated in the discussion of the "inner" solution, the effect of 
changing initial wave amplitude Lo from 0.1 to 0.01 is indistinguishable on 
the scale of Fig. 4.1, except for the level of L between £ = 0 and 2, where the 
drop appears. 

One sees that when Q passes through zero again, at about £ = 7 , 
another sudden drop occurs, with a reappearance of rotating-wave amplitude 
L. The reason for this event is that, following the first drop (at about C, = 2), 
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Q is negative, and the rotating wave is damped exponentially to a very small 
amplitude; but then the surge effect provides enough re-amplification 
(between £ about 3 to 7), so that a non-zero resonance occurs when Q 
passes down through the axis again. This second "inner" type of wave is 
weaker than the first because the corresponding velocity drop is weaker, as 
Eq. (4.40) implies. 

It should be emphasized that the strength of the surge oscillation shown 
in Fig. 4.2 is due to the resonance. Fig. 4.3 shows how the transient would 
develop if there were no resonance effect at all. (In the calculations, L was 
held equal to zero.) A very weak surge is indicated, as flow coefficient 
follows the changing throttle coefficient. 

Effect of throttle closure rate . Figs. 4.4 and 4.5 show the 
corresponding results if the transient begins at a much higher flow 
coefficient, and the throttle closure is much faster. Qualitatively, the results 
are the same for the surge effect, and the resonant drop is quantitatively the 
same, as expected, because it depends only on conditions existing just when 
resonance occurs. 

Figs. 4.6 and 4.7 show results when the throttle is abruptly closed at £= 
0, to a level calling for a change of flow coefficient from 1.1 to -1.1, which is 
below the resonant value. Fig. 4.6 shows that, again, a surge develops after 
the sudden transition, but there is no second resonance, because the surge 
does not penetrate into the amplifying range of Q between 0 and 1. Fig. 4.7 
shows the contrasting result if the rotating wave L is suppressed. In this 
case, the surge is stronger. In effect, resonance gives a weaker surge because 
the sudden drop happens to place Q close to the level demanded by the new 
throttle setting anyway. 

Fig. 4.8 shows what happens when the throttle is abruptly closed again 
from 1.1, but only to 0. As might be expected, after the first sudden drop, 
the flow coefficient hunts around the resonant value, with an indefinite 
succession of weak resonant transitions. 

Effect of characteristic slope . Next, the effect of characteristic slope 
(a) is considered. One expects that positive values will lead to amplification 
of surge, and perhaps of rotating waves as well; negative values typically 
damp surge, and they certainly suppress rotating stall in axial machines. So 
far, only the neutral situation (o = 0) has been discussed. 

First, a small positive (destabilizing) slope of <T = 0.2 is considered, 
with a slow throttle closure. The result, shown in Fig. 4.9, may be compared 
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with Fig. 4.2, for which o = 0. The result is similar, except that the surge 
induced be the resonance is amplified, as expected. Also, there seems to be 
no secondary resonance, presumably because the rotating-wave damping is 
so strong. 

Fig. 4.10 then shows the result when the characteristic slope is 
negative. The surge motion is indeed damped as one expects, but secondary 
resonance events of some strength appear. If the slope is more strongly 
negative, as in Fig. 4.1 1, this tendency is accentuated; a succession of 
resonant transitions occur which are about as strong as the first one, and they 
occur in rapid succession. The flow coefficient shows a saw-tooth pattern 
with time; a sort of "stick-slip" oscillation about the throttle value. The 
pressure-rise coefficient (Z) drifts toward negative values, which is opposite 
to the trend expected when the throttle closes and the characteristic slope is 
negative. Fig. 4.1 1 may be contrasted with Fig. 4.12 which is the same 
except that rotating waves are suppressed. In the latter case, surge-like waves 
are weak, and the pressure rise ultimately begins to rise in response to the 
throttle. Clearly, resonance is predicted to produce a dramatically different 
transient behavior, even when on the "stable" slope of the characteristic. 

The successive surge transients shown in Fig. 4.1 1 continue until the 
peaks of Q become negative; then the rotating waves are no longer 
amplified. Fig. 4.13 illustrates this effect for a higher throttle closure rate. 
After about C, = 5, the resonance no longer operates, because Q stays below 
1 . Subsequently, Q follows the throttle more closely, and Z steadily 
increases as the characteristic slope requires. 

Finally, a still steeper "stable" slope of a = 1 gives the results of Fig. 
4.14, in which the tendencies just described are accentuated. 

Resonant capture of surg e. If a surge transient occurs about a throttle 
setting which is above the range of amplification, and the surge happens to 
carry Q into the amplification range, the sigularity "captures" the 
oscillation. Figs 4.15, 4.16, and 4.17 show this effect, for a strongly stable 
a = 1. The throttle is set for 1.1, 1.4, and 2.4, respectively. In each case, 
begins at 2.0. Fig. 4.15 shows successive capture, which presumably recurs 
indefinitely; Fig. 4.16 shows capture, but insufficient time in the 
amplification range for recurrence; Fig. 4.17 shows a case with the throttle 
setting so high that the amplification range is not entered at all, and there is 
no capture. 
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Resonance for Higher Harmonics 

Fig. 1.2 was the basis for the foregoing analysis, and it concerned the 
first circumferential harmonic for a traveling wave. The singularity was 
associated with the lower branch of the wave-speed (f) versus flow 
coefficient (Q) function, and f would be quite stationary at zero, near the 
singularity, for the particular case chosen for study. Actually, as explained 
in [1], the third or fourth harmonics are more technically interesting; they 
would be found at higher flow coefficient according to both theory and 
experiment. 

The shape of the amplification curve for the third harmonic, for 
example, is different from that shown in Fig. 1.2; Fig. 2.2 suggests that the 
function f(Q) has a cusp at the singular point, and that amplification occurs 
for the upper branch of that function. There is virtually no damping on that 
branch; damping is associated with a slightly lower wave speed. 
Unfortunately, there is no easy way to evaluate singularity strengths in this 
case, because the simplification of f = 0 is not available for any reasonable set 
of parameters. Nevertheless, one may speculate that equations similar to 
Eqs. (4.33) through (4.35) would still apply, after certain changes are made: 
For negative Q, amplification is positive; therefore, the sign of the L/Q term 
of Eq. (4.34) is changed. It is assumed that the last term of Eq. (4.33) 
remains as written, because it depends on the square of Q. Further, it is 
imagined that both those terms apply only when Q is between 0 and -4; 
otherwise, they are both taken to be zero. 

The solution of the revised equations is quite different from the 
previous one. The "inner" solution 

Q = -|C; Z = 0 ; L=|C 2 (4.47) 

applies even with no initial value of L. The maximum L is 8, whereas it was 
only 0.5 in the previous case. Also, because there is no damping, wave 
amplitude remains high even though Q may pass below the range of 
amplification; previously, L would return to zero. Fig. 4.18 shows a 
numerical calculation, on the "inner" scale, for the response to a sudden drop 
of throttle flow. 

Fig. 4.19 follows the response on the long-term scale (Fig. 4.6 shows a 
comparable previous case.) Following the initial transient, a surge develops 
as before, but now the wave amplitude does not return to zero; in fact, the 
next time the surge comes to Q = -1, a new increment of L appears. As one 
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would expect, if the characteristic slope is even slightly positive 
(destabilizing for surge), this effect is accentuated (Fig. 4.20), and with each 
surge cycle, the traveling wave would be kicked to a considerably higher 
level. If the slope were negative, no such kick would occur, and L would stay 
at the level of 8. 
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5. TRANSIENTS WITH MULTIPLE ROTATING MODES 

Away from resonant conditions, the possibility of more than one wave 
type should be considered, because two families of rotating waves (fast and 
slow) are in principle allowed. Furthermore, many different slow-wave 
modes could possibly coexist. That is, instead of Eqs. (3.10) and (3.1 1), one 
would begin with 

g = W [Aj 1} ) cos (9 - fO)^) + A (2) ) cos (0 - f( 2 )^) + ...] (5.1) 


P = P (r,£) + iP( 1 )cos(0 - fU)£) + 2P (1) sin(0 - f(U£) 

+ 1 P(2)cos (0 - P 2 )^) + 2P^2) sin(0 - tf 2 )^) + . . . (5.2) 

It seems clear that the Galerkin procedure can be applied only to angle 
variation, not time; the various modes are not orthogonal, and therefore 
averages over time cannot be defined. If appropriate expressions based on 
Eqs. (5.1) and (5.2) are introduced into Eq. (2.2), and then angle averages 
and the cos0 and sin0 moments are derived, a very complex set of unsteady 
equations result. They were not solved; however, they imply that one should 
expect surge-like beat pulsations depending on frequency differences, and 
one should expect some exchange effects between modes. These things would 
depend on products of amplitudes, of course, and would therefore be weak 
unless both component waves were quite strong. 

It seems generally unlikely that two or more circumferential modes 
would be present in substantial strength at the same time. The various 
possible modes are subject to separate and different rates of amplification, 
often with resonance for some particular mode. Thus, it seems that one mode 
will almost always be singled out by the amplification process, and will 
predominate over the others. 

For these reasons, the problem of multiple rotating modes is not 
pursued further in this report. 
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6. CONCLUSIONS 


(1) The previous small-perturbation analysis of circumferential 
waves, assuming inviscid flow, does not depend on diffuser width; however, 
it does allow for radial convergence of the diffuser. Convergence confers 
stability in that the critical flow coefficient for the appearance of waves is 
reduced. The observation that narrow diffusers are more stable than wide 
ones is argued to be the result of boundary-layer displacement effect tending 
to make narrow diffusers relatively more convergent than wide ones. 

(2) The Moore-Greitzer analysis for finite-strength transients 
containing both surge and rotating stall in axial machines can be adapted for 
the centrifugal compression system. In the axial case, the motivation was to 
explain "hung stall". In the centrifugal case, the motivation is to predict 
events when flow is throttled through resonance. As previously, a Galerkin 
procedure is adopted, and a cubic characteristic shape is assumed. 

For a transient, it is assumed that wave speed is constant, that only 
squares of wave amplitude need be considered, and that flow in the diffhser, 
though time-dependent, is quasi-steady. These assumptions permit nonlinear 
expressions to be derived for pressure change and for a new, second-order 
swirl velocity. Except for this swirl, the diffuser velocity field can be 
carried over from the previous linear analysis. 

The nonlinear equations confirm that negative values at the diffuser 
entrance of the coefficient of the pressure wave make the circumferential 
wave amplitude grow with time. 

(3) The singularities of interest are simple poles in general. The pole 
strengths can be found analytically for a particular case of the fundamental 
mode and radius ratio of about 2; for that case, propagation speed happens to 
be zero. 

Using the nonlinear transient equations, it is shown that throttling 
through a resonant value of flow coefficient cannot be smooth; flow 
coefficient must suddenly drop through that value in a surge-like transient, 
accompanied, of course, by a transient rotating wave. As in the axial 
analysis, it is necessary to postulate some pre-existing low initial level of 
circumferential disturbance, though the particular value chosen turns out not 
to be important. 

A simple example shows this sudden drop to be an "inner" solution of 
the transient equations, while an "outer" solution on a longer time scale 
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governs the surge response to the initial excursion of flow coefficient. 
Typically, the long-term surge is much stronger than would be expected 
from throttle change alone. In fact, if throttle is set above the critical flow 
coefficient, the singularity may "capture" the flow, and a surge would occur 
purely as result of circumferential-wave resonance. 

Numerical results are shown for various parametric choices relating to 
the throttle schedule and the characteristic slope. Especially interesting is the 
result of ramping through resonance when the characteristic slope happens to 
be negative; a rapid succession of surge excursions occur in a saw-tooth 
pattern. 

(4) For higher circumferential modes, the simple pole type of 
singularity occupies a very narrow range of flow coefficient, and strong 
amplification seems to dominate for a range of flow coefficient beginning 
just below the critical. No analytical solution is available for such a case. 
However, it is suggested that the form of the result will involve strong surge 
oscillations, accompanied by rotating waves which do not die out, but rather 
accrete strength every time the surge passes through the resonant level. 

(5) In principle, a number of circumferential modes should be 
considered simultaneously, because they all might be candidates to participate 
in a transient. It is shown how such a general Galerkin analysis would 
proceed, and it is expected that surge-like beat oscillations and mode- 
exchange effects could arise. However, it is argued that modes will probably 
appear separately because they are separately amplified, especially if 
resonance is involved. Thus, multiple circumferential modes are not treated 
further. 

(6) To construct exact solutions of Eqs. (3.42)-(3.45), one would first 
establish a file corresponding to Figure 2.2, for example, including velocity 
and pressure profiles. That information would constitute the quasi-steady 
solution to be used as O varies. Eq. (3.45) would have been satisfied in 
forming that solution. Next, one would solve the remaining equations by 
marching in time; at each time step, the integral in Eq. (3.42) would be 
evaluated. Every time that O descends through a resonant point, the 
appropriate "inner" transition would be applied. 
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Fig. 1.1 Sketch of centrifugal compression system, schematically comprising 
entrance, rotor, diffuser, plenum, and throttled discharge. 
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Figure 1.2 Amplification and wave speed for 1st mode (b=l), showing effect of 
radius ratio r p . Negative values of PRj signify amplification. 
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Figure 2. 1 Sketch contrasting the boundary-layer displacement effect 
for a narrow and a wide diffuser. 
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Figure 4.4 Effect of rapid throttle closure; otherwise same as Figure 4.2. 
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Figure 4.5 Same as previous Figure, except that resonance is suppressed by 
holding L= 0. 
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Figure 4. 19 Development of surge for amplification typical of higher harmonics, with 
zero (neutral) characteristic slope. 
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Figure 4.20 Development of surge disturbance for amplification typical of higher 
harmonics, for a slightly negative (stable) characteristic slope. 



APPENDIX A 


In [ 1 ], owing to length limitations, the solution of the small- 
perturbation equations was presented only for parallel diffuser walls (n = 0) 
and the first circumferential mode (b = 1). For the record, the general 
solution for any n and any b is described in this Appendix. The 
correspondence with equations in [1] will be noted, and, where possible, the 
same notation will be used, with the caret ( A ) applied to specify quantities 
arising in that linear analysis. However, in this report, O rather than V is 
used to denote mean flow coefficient. Certain new quantities j, x, y, p, and v 
appear only in this Appendix; they are defined in the text below. 

The applicable inviscid equations for flow in the diffuser are recorded 
in Eqs. (3.17)-(3.19). Small perturbations are assumed, in terms consistent 
with Eqs. (3.20) and (3.21), and with Eqs. (6)-(8) of [1]: 


v O 
U" r 


A e ibe 


^ = Weib9 ; d P = Pe i b0 


(A.l) 


The three Euler equations for small disturbances then become, instead of 
Eqs. (9)-(ll) of [1], 

Vr" [rA r - (1 - n) A] + ib (1 - fr2)A - 2W + = 0 (A.2) 


Vr n (rW r + W) + ib (1 - fr2) W + ir P = 0 


(A. 3) 


r A r + (1 -n)A + ibW = 0 (A.4) 

The boundary conditions corresponding to Eq. (12)-(14) of [1] are now, at 
the diffuser entrance (r = 1), 

A(l) = 1 ; W (1) = 0 (A. 5) 

and at the diffuser exit (r = r F ), 

P (rF) = 0 (A. 6) 

The foregoing equations reduce to those in [1] if n and b are set equal to 0 and 
1 , respectively. 
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The solution of this system of equations for P, A, and W has the 
following form (only P was given in Eq. (15) of [1]): 

b P = [Ki Pi . K 2 P 2 - PPl] + i [K 2 Pi + Ki P 2 - PP 2 ] (A.7) 

A = r [Ki (yij 3 - yzjl) - ^2 (yij4 - yzj2>] + (p + 1 - n) ryi 

+ (|4 + 1) r j 2 + ir [K 2 (yij 3 - Y 2 jl) + Kl (yij4 - y 2 J 2 )] (A.8) 

bW = Ki [(p + 1) yij4 + (P + 1 - n) yzJ 2 ] + K 2 [(p + 1) yij3 + 

(p + 1 - n) yzji] + i { - Kl [(p + 1) yij3 + (p + 1 - n) yzji] + 

K 2 [(ft + 1) y U4 + (p + 1 - n) yzj 2 ] + r (y 2 - yi) } (A.9) 

where, corresponding to Eqs.(16)-(18) of [1], the following definitions 
apply: 

Pi = ^ [cos z - b 2 r n (y ij3 - yzjl)] - (1 - fr 2 ) [(p + 1 ) yij4 


+ (p + 1 - n) y 2 j 2 ] 

(A. 10) 

P2 = ^ [sin z - b 2 ^ (yij4 - y2J2)] 


+ (1 -fr 2 ) [(p + l)yij3 + (P+ 1 - n) yzii] 

(A. 11) 

PPi = ^ r n b 2 [(p + I)y 2 + (P+ 1 -n)yi] 

(A. 12) 

pp 2 = 1 ) 2(1 - fr 2 ) (y 2 - yi) 

(A. 13) 

Two parameters involve n and b: 


1 n ~lu7 n2 

P=- 1+2 + ^/ b2 + 4 " 

(A. 14) 

v = n-3-2p 

(A. 15) 


and certain functions have the definitions which follow. 


I 


56 


11 1 mu 1 1 l 1.1 i 111 iin 11 1 1 11 min 


z = (/nr + ^fr 2 ) ^ if n = 0 


(A. 17) 


yi = 

r v+p 
" v+ 1 

(A. 18) 

y 2 = 

rP-i 

" v+ 1 

(A. 19) 

It may be noted that z is the same as $ defined in Eq. (4.8) when f = 0. 

The following integrals require numerical evaluation, in which x 
denotes the variable of integration from 1 to r: 

jl = 

r 

r A 

J x'(l I+n ) cos z(x) dx 

1 

(A. 20) 

J2 = 

r 

f A 

J x*(l I+n ) sin z(x) dx 
1 

(A. 21) 

j3 = 

r 

J x ‘( 3n " 1 + v V 2 cos z(x) dx 
1 

(A.22) 

j4 = 

r 

J x *(3n - 1 + v)/2 s j n z ( x ) dx 

(A.23) 


1 


Inspection shows that the solution just described automatically satisfies 
the boundary conditions (Eq. (A.5) at the diffuser entrance. Those 
conditions require modification [1] if the rotor blades are not radial at the 

A A 

rotor exit, of course. The constants Ki and K 2 may be chosen to satisfy the 
boundary condition at the diffuser exit, Eq. (A.6), with these results: 



(A. 24) 


A AA AA A A 

Ki = [(Pi PP1 + P 2 PP2)/(Pl 2 + P2 2 )] r = rp 
Kl = [(Pi PP2 - P2 PPl) /(Pi 2 + P2 2 )] r = rp (A. 25) 

The latter boundary condition is changed if there is any recovery of dynamic 
head at the diffuser exit, as explained in [1]. 
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NOTATION 


The notation in this Report is consistent with that in [1] and [4], because 
the analysis is, in effect, an amalgamation of those previous two papers. All 
lengths are in ratio to the radius of the diffuser inlet (rotor discharge). 

A amplitude of cosine part of circumferential wave of radial 
velocity disturbance, Eqs. (3.10), (3.20) 

Af effective flow area 

AM imaginary part of disturbance A in [1], Eq. (3.35) 

AR real part of disturbance A in [1], Eq. (3.35) 
a reciprocal time-lag parameter of blade passage, Eq. (3.4) 

B amplitude of sine part of velocity disturbance, Eq. (3.20) 

B g inertia parameter for surge ("Greitzer's B") 
b mode number for circumferentially-traveling wave 

b* pole strength for pressure disturbance, Eq. (4.20) 
c* pole strength for integral, Eq. (4.20) 

Ft throttle characteristic function, Eq. (3.8) 

f coefficient of propagation speed of circumferential disturbance, in 
laboratory frame, Eq. (3.12) 

G function of exit radius and flow coefficient defined in Eq. (4.9) 
g disturbance of axial flow coefficient, Eq. (3.5) 

H height parameter of cubic axisymmetric characteristic. Fig. 3.1 
I integral expression defined in Eq. (4.17) 

J square of disturbance wave amplitude at diffuser entrance, Eq. (3.39) 
Ko ititial throttle setting, Eq. (4.46) 
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Ki rate of throttle closure, Eq. (4.46) 
k constant in "inner" solution, Eq. (4.41) 

L amplitude of circumferential disturbance of radial velocity, Eq. (4.32) 
Lo initial value of L 

/ diffuser width, Eq. (2.1) and Fig. 2.1; generally a function of radius 
/j length of inlet section before rotor. Fig. 1 of [4] 
l c effective flow length of system up to rotor exit, Eq. (3.4) 

M function of exit radius and flow coefficient defined in Eq. (4. 10) 

m duct-flow parameter of rotor entrance, Eq. (3.6) 

n exponent describing diffuser convergence, Eq. (2.1) 

p static pressure in diffuser 

Pf plenum pressure to which diffuser discharges 
Pt upstream total pressure 

P pressure coefficient, Eq. (3.7) 

P angle-averaged pressure coefficient, Eq. (3.7) 

Pqs quasi-steady pressure coefficient 
d P wave-like part of pressure coefficient, Eq. (3.7) 
iP amplitude of cosine part of pressure wave, Eq. (3.11) 

2 P amplitude of sine part of pressure wave, Eq. (3.1 1) 

PM imaginary part of pressure-wave amplitude from [1] 
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PR real part of pressure- wave amplitude from [1] 

Q nonlinear disturbance of flow coefficient, Eq. (4.32) 

q nonlinear disturbance of flow coefficient, Eq. (4.24) 

r radial coordinate (relative to radius of diffuser inlet) 
rp radius at diffuser exit; radius ratio of diffuser, Fig. 1.1 
S shut-off parameter of cubic axisymmetric characteristic, Fig. 3.1 
SL local slope of compressor characteristic, (Eq. (4.21) 
t time 

U wheel speed at rotor exit (diffuser inlet) 

v radial velocity in diffuser 

W width parameter of cubic axisymmetric characteristic. Fig. 3.1 
w circumferential velocity in diffuser 

X second-order swirl, Eq. (3.21) 

Y potential for wave-like velocity disturbance, Eq. (3.5) 

Z nonlinear disturbance of pressure-rise coefficient, Eq. (4.32) 

z nonlinear disturbance of pressure-rise coefficient, Eq. (4.25) 

a coefficient in analysis for f = 0, Eqs. (4.2)-(4.4) 

P coefficient in analysis for f = 0, Eqs. (4.2)-(4.4) 

y coefficient in analysis for f = 0, Eqs. (4.2)-(4.4) 

e initial small disturbance of flow coefficient, Eq. (4.32) 

X strength of wave disturbance, (Eq. (4.26) 
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d> flow coefficient based on velocity at rotor exit and wheel speed there 
*F actual pressure-rise coefficient, Eq. (3.3) 

\j/ c axisymmetric characteristic function, Eq. (3.13) 

\|/co height parameter of axisymmetric characteristic function, Eq. (3.13) 
0 circumferential angle coordinate, laboratory frame 

0 angle coordinate in frame of traveling wave, Eq. (3.12) 

$ function of exit radius and flow coefficient, Eq. (4.8) 
p density 

a redefined slope of characteristic, Eq. (4.28) 

x blade-passage lag coefficient 

£ time coordinate; rotor travels one radian in unit time 
£ redefined time coordinate, Eq. (4.27) 

"inner" time variable, £/e 

Superscripts 
* at resonance 

A obtained from the linear analysis of [ 1 ] 

0) identifies traveling mode in Eqs. (5.1), (5.2) 

(2) identifies traveling mode in Eqs. (5.1), (5.2) 

primes denote differentiation in Eqs. (4.14) and (4.15) 
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Subscripts 


o at rotor entrance 

j at diffuser entrance 

F at diffuser exit 

T pertaining to throttle (except pj, which means ambient toal pressure) 
partial differentiation denoted by subscripts 


63 



REFERENCES 


1 . Moore, F.K. "Weak Rotating Flow Disturbances, in a Centrifugal 
Compressor With a Vaneless Diffuser," ASME J. Turbom achinery, 111, 
October 1989, pp. 442-449. 

2. Moore, F.K. "A Theory of Rotating Stall of Multistage Compressors: 
Part I— Small Disturbances," ASME J. Eng, for Power , 106, April 1984, pp. 
313-321. 

3. van den Braembussche, R., "Rotating Stall in Vaneless Diffusers of 
Centrifugal Compressors," Tech. Note 145, von Karman Institute for Fluid 
Dynamics, Rhode Ste Genese, Belgium, June 1982. 

4. Moore, F.K., and Greitzer, E.M., "A theory of Post-Stall Transients on 
Axial Compression Systems: Part I— Development of Equations," ASME L 
Eng, for Gas Turb. and Power , 108, January 1986, pp 68-76. 



Report Documentation Page 

Space Administration 

1. Report No" 2. Government Accession No. 

NASA CR- 1 87049 

4. Title and Subtitle 

Theory of Finite Disturbances in a Centrifugal Compression System with a 
Vaneless Radial Diffuser 


3. Recipient's Catalog No. 


5. Report Date 
December 1990 

6. Performing Organization Code 


7. Author(s) 

F.K. Moore 


9. Performing Organization Name and Address 

Cornell University 
Ithaca, New York 14853 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

15. Supplementary Notes 


8. Performing Organization Report No. 


10. Work Unit No. 

535-05-10 

11. Contract or Grant No. 

NAG3-349 

13. Type Qf Report and Period Covered 
Contractor Report 

Final 

14. Sponsoring Agency Code 


Project Manager, Daniel H. Buffum, Propulsion Systems Division, NASA Lewis Research Center. 


16. Abstract 


A previous small-perturbation analysis of circumferential waves in circumferential compression systems, assuming 
inviscid flow, is shown to be consistent with observations that narrow diffusers are more stable than wide ones, 
when boundary-layer displacement effect is included. The Moore-Greitzer analysis for finite-strength transients 
containing both surge and rotating stall in axial machines is adapted for a centrifugal compression system. Under 
certain assumptions, and except for a new, second-order swirl, the diffuser velocity field, including resonant 
singularities, can be carried over from the previous inviscid, linear analysis. Nonlinear transient equations are 
derived and applied in a simple example to show that throttling through a resonant value of flow coefficient must 
occur in a sudden surge-like drop, accompanied by a transient rotating wave. This inner solution is superseded 
by an “outer” surge response on a longer time scale. Surge may occur purely as result of circumferential-wave 
resonance. Numerical results are shown for various parametric choices relating to throttle schedule and the 
characteristic slope. For higher circumferential modes, a somewhat different picture of the amplification process 
seems appropriate, and suggestions are made about the effects of resonance in such cases. The problem of 
considering a number of circumferential modes simultaneously is briefly discussed. 


17. Key Words (Suggested by Author(s)) 

Compressor; Stability; Centrifugal; Stall; Surge; 
Resonance; Diffuser 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 07 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of pages 
66 



NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 






